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We investigate whether or not an Einstein Static universe is a solution to the cosmological equa- 
tions in f{R) gravity. It is found that only one class of f{R) theories admits an Einstein Static 
model, and that this class is neutrally stable with respect to vector and tensor perturbations for 
all equations of state on all scales. Scalar perturbations are only stable on all scales if the matter 
fluid equation of state satisfies c^ > ""*" ~ 0.21. This result is remarkably similar to the GR case, 



o 

O 



(N 



I . 



> 

00 
(N 

vn 

en 

^' 

o 

oo 

O 



X 



PACS numbers: 04.50.+h, 04.25.Nx 



I. INTRODUCTION 

Although modifications of Einstein's theory of gravity 
were already proposed in the early years after the pub- 
lication of General Relativity (GR), a detailed investiga- 
tion of cosmological models within this framework only 
got underway a few years ago. Such models became pop- 
ular in the 1980's because it was shown that they natu- 
rally admit a phase of accelerated expansion which could 
be associated with an early universe inflationary phase 
[l|. The fact that the phenomenology of Dark Energy 
requires the presence of a similar phase (although only 
a late time one) has recently revived interest in these 
models. In particular, the idea that Dark Energy may 
have a geometrical origin, i.e., a that there is a connec- 
tion between Dark Energy and a non-standard behavior 
of gravitation on cosmological scales is currently a very 
hot topic of research. 

In the last few years cosmologies based on fourth order 
gravity have been widely studied using several different, 
but complementary techniques. A number of extremely 
interesting results have been found 0, 0, \^ which sug- 
gest that it might be possible that these models pro- 
vide a viable alternative to Dark Energy [^, [g, 0, Q. 
Among these approaches, a particularly interesting one 
is based on the analysis of the phase space of these mod- 
els 0, [13, Uil [H) [3' providing a systematic method 
for obtaining exact solutions together with their stability 
and a general idea of the qualitative behavior of these 
cosmological models. 

Of particular interest is the fact that there are classes of 
f{R) theories which admit a Friedmann transient matter- 
dominated decelerated expansion phase, followed by one 
with an accelerated expansion rate [0, Q . The first phase 
provides a setting during which structure formation can 
take place and this is followed by a smooth transition 
to a Dark Energy like era which drives the cosmological 
acceleration. It would therefore be of great interest if 
orbits could be found in the phase-space of f{R) models 
that connect such a Friedmann matter dominated phase 
to an accelerating phase via an Einstein Static solution in 
a way which is indistinguishable (at least at the level of 



the background dynamics) to what occurs in the ACDM 
cosmology of GR [ij. If a similar evolution exists in 
an /(i?) model, we would expect the associated Einstein 
Static solution to be a saddle point (as it is in GR) and 
consequently unstable. 

The issue of stability of the Einstein Static model may 
at first seem a somewhat trivial matter. Indeed, already 
in the 1930s, Eddington 15] showed that in GR such 
models are unstable with respect to homogeneous and 
isotropic perturbations, and ever since then the Einstein 
Static model has been considered unstable to gravita- 
tional collapse or expansion. Indeed, it is exactly this 
feature that allows the transition between a decelerated 
expansion era to one which is accelerating in the ACDM 
cosmology. However, later work first by Harrison [l6| 
and then by Gibbons [17] suggests that this issue requires 
further investigation, both in the context of GR and ex- 
tended theories of gravity. In Harrison's classic paper, 
he investigated the normal- modes of vibration of the uni- 
verse and showed that all physical inhomogeneous modes 
are oscillatory in a radiation-filled Einstein Static model, 
while Gibbons showed that an Einstein Static model filled 
with a perfect fluid was stable with respect to conformal 
metric perturbations, provided that the sound speed sat- 
isfies c^ > i The reason for this "non-Newtonian" sta- 
bility stems from the fact that in GR, the Einstein Static 
universe is spatially closed therefore has a maximum scale 
associated with it, which is greater than the largest physi- 
cal scale of the perturbations. Since the Jean's length is a 
significant fraction of this maximum scale, perturbations 
in the fluid oscillate, rather than grow, leading to con- 
clusion that in GR at least, the Einstein Static solution 
is marginally stable with respect to such perturbations. 
This result was also found more recently by Barrow et. 
al. [18], when linear inhomogeneous and anisotropic per- 
turbations of the Einstein Static model were investigated. 

Another reason why Einstein Static solutions and their 
stability are both interesting and important is due to 
the fact that a number of authors have conjectured that 
the universe might have started out in an asymptotically 
Einstein Static state, providing a natural beginning to 
the inflationary expansion phase of the Universe [101 . 



The issue of stability also plays an important role in 
the viab ility of cosmological solutions in f{R) theories 
[13, lU, [11 H- Recently [H it was argued that the 
sign of the second derivative f^RR = d^ f /dB? determines 
whether the theory approaches the general relativistic 
limit at high curvatures, and it was shown that for /,_r_r, > 
the models are, in fact, stable. The stability of the de 
Sitter solution in f{R) gravity has also been extensively 
analyzed in the literature [2a |. 

Motivated by the above discussion, the goal of this pa- 
per is to first determine under what conditions an Ein- 
stein Static solution exists in a general f{R) gravity the- 
ory and investigate the stability of such solutions with 
respect to general inhomogeneous and isotropic pertur- 
bations. 



II. EINSTEIN-STATIC UNIVERSES IN 
FOURTH ORDER GRAVITY 

In a completely general context, a fourth order theory 
of gravity is obtained by adding terms involving RabR"^'' 
and Rabcd.R"''"^'^ to the standard Einstein Hilbert action. 
However, it is now well known that if we use the Gauss 
Bonnet theorem we can neglect the Rahcd.R"'^^'^ term [g^ . 
Furthermore, if we take into account the high symmetry 
of the FRW metric, the Lagrangian can be further sim- 
plified. Specifically the variation of the term R^bR"^^ can 
always be rewritten in terms of the variation of R^ I27l| . 
It follows that the " effective" fourth-order Lagrangian in 
FRW cosmology only contains powers of R and we can, 
with out loss of generality, write the action as 



A = 



d^x^ 



-g[!{K)-2K + C,n] 



where Cm represents the matter contribution and A is 
the usual cosmological constant. The corresponding gen- 
eralization of Einstein's equations are 



r Jl, „ ^^ _ rpT _ ±ab 
•^ab + yab-p — J^ ab — "FT 



T, 



ab 1 



(2) 



where /' = f{R),R. Here T^ is the total effective energy 
momentum tensor, T^'jJ is the energy momentum tensor 
for standard matter and 



rpR. 

^ab 



1 

J 



1 



UI-Rn^fU9'a9t-g'''9ab) 



(3) 



is the energy momentum tensor containing all the non- 
GR curvature contributions, known as the ^^ curvature 
fluid' . 

The metric for a closed Einstein Static universe is given 
by 



ds^ = -di^ + al 



dr"^ 



l~r^ 



r'^dn^ 



(4) 



Here ap is a constant and both the Ricci scalar R and the 
curvature of the 3-space R are equal and can be written 



as 



R — R — —77 



a, 



(5) 







It is easily seen that an Einstein Static universe is expan- 
sion, shear and rotation free, i.e. 

6 = 0, aab - 0, c^ab = . (6) 

The total energy momentum tensor is given by 



T 1 

T 1 






(7) 



where pm and Pm is the energy density and pressure of 
standard matter. 

Now the (t) component of Einstein's equations gives 



R 

"2 



P 



A 



while the trace of Einstein's equations is 

A 



R 



/' 



3P^ 



(8) 



(9) 



Furthermore, from the Raychaudhuri equation we have 
Rabu'^u^ = 0, from which we obtain 



(1) Using 



p' + 2>p' - 

and pO|) we find 

i? = 3 [p^ 



A 



p 



(10) 



(11) 



Now using the components of total energy momentum 
tensor ([7]), and assuming that the standard matter is a 
barotropic perfect fluid with equation of state and sound- 
speed parameters w and cj. defined by 



p" = wp"", -1 < w < 1 
equation (fTTj) becomes 

Rf = 3(1 + w)p" 



w , 



(12) 



(13) 



The LHS of the above equation cannot be zero for this 
spacetime, because this would imply either oq — > oo or 
the function f{R) — constant. Consequently, we cannot 
have an Einstein Static universe without matter or with 
w = —1. 

It follows from p^ above, together with the trace 
equation ^ and ([7]) that for an Einstein Static universe 
to exist, the function f{R) needs to satisfy the following 
differential equation: 



Rf'-T;{^ + w)f + 3A{l + w)^0 



(14) 



In order to obtain the GR limit we set f = R and f — 1. 
In this case we obtain 



R 



6 _ 6A(l + u;) 



(15) 



It is easy to see that in GR, once we fix oq and A, both 
w and p™ get fixed via equations p^ and (fT3|) . The 
general solution for (I14p is given by 



/(i?) = 2A + /Ci?i(i+'"' 



(16) 



where /C is the constant of integration. It follows that 
the action for the Einstein Static universe ^ becomes 



A^ 



/Ci?5(i+"')+£, 



(17) 



We can then easily see that in fourth order gravity, the 
existence of an Einstein Static universe implies there is 
no explicit dependence on the cosmological constant in 
the action. However, we know that in GR an Einstein 
Static universe exists for the action 



A^ 



^[R-2A + C„ 



(18) 



with the Ricci scalar given by ([15)) . As noted earlier, 
the basic difference between GR and fourth order the- 
ories of gravity, in case of Einstein Static Universe is, 
for GR fixing w and A would fix p™ and Uq simultane- 
ously. Whereas in f{R) gravity fixing the former quan- 
tities would fix the 'function' / and not the values of 
p"* and R. Hence we can in principle have an Einstein 
Static Universe with any non-zero positive p™ for a given 
barotropic index and cosmological constant and this is a 
direct consequence of the extra degrees of freedom for 
higher order gravity. However if we wish to impose the 
condition that for the value of the Ricci scalar given by 
(fT5|) . the function f{R) also has the same value (so that 
we have a GR limit), that would then fix the constant of 
integration K. as 



/C 



4A 



{l + 3w) 



{l + 3w) 



-I i(i+«') 



6A(l + 'u;) 



(19) 



However, if we switch off the cosmological constant term 
in ([T]) from the very beginning, it is possible to get an 
Einstein Static universe for any value of /C. 

It is extremely interesting to note that just the back- 
ground symmetry of the Einstein Static universe fixes the 
form of allowed actions in the fourth order gravity the- 
ory. It follows that, based on the prescription leading to 
equation (14), many commonly used forms of fourth or- 
der gravity, such as /(i?) = R + aR" or f{R) = R + a/R 
do not admit an have an Einstein Static solution. 



III. STABILITY ANALYSIS 

In the previous section we described the complete dy- 
namical equations describing an Einstein Static universe 



in fourth order gravity and discovered that there is a 
unique class of f{R) theories that admits such a solu- 
tion. In what follows, we investigate the stability of this 
model with respect to generic linear inhomogeneous and 
anisotropic perturbations by linearizing the most gen- 
eral propagation and constraint equations for this f{R) 
theory around the Einstein Static background. This is 
done using the 1-1-3 covariant approach to perturbations 
[28l . [29l . [so . [31I [33 . [33 | , where quantities that vanish in 
the background spacetime are considered to be first or- 
der and are automatically gauge-invariant by virtue of 
the Stewart and Walker lemma 3j|. In order to write 



down the set of linear equations we first need to choose 
a physically motivated frame u°. The natural choice is 
the one which is tangent to the matter flow lines. 

The projection tensor into the tangent 3-spaces orthog- 
onal to this flow vector is 



hab = gab + UaUb 



(20) 



Also, the totally projected covariant derivative operator 
Va, orthogonal to u° is defined by 



Va = K^b 



(21) 



The effective total density and pressure in the presence 
of perturbations is then 



T 1 



T ^ 



p 



\{Rf-f) + I"^^R 



liRf'~f)^lf"y'R + f"R 



(22) 



(23) 



while the anisotropic part of effective energy momentum 
tensor is given by 



H^-1 

"' ' r 

and the energy flux is 



/"V<aV6>i? 



T 1 



r'RVaR+f"yaR~^f"VaR 



(24) 



(25) 



Here angle brackets denote the projected trace-free part 
(see [33 for details). 

In terms of these effective thermodynamical quantities, 
the perturbed cosmological equations are 



e^ 



a = gBa , 



^ A' 

p^ + j 



1^ 



(26) 



(27) 



In addition to these, energy conservation for standard 
matter is given by 



-ep"\i + w) 



(28) 



while the momentum conservation gives the following re- 
lation connecting the acceleration to the matter density: 



A" =u'' = — 



w + 1 p™ 



(29) 



The expansion propagation (or the generalized Ray- 
chaudhuri equation) is 

e = -ie^ + v'^A. -\{p^ + 3p1 + J ■ (30) 

The shear propagation equation can be written as 
2 1 ~ 

CTab = -^Oo-afc - ^ah + -^nafc + V(Q^h) , (31) 

where the tensor Eat is the "electric" part of the Weyl 
tensor. The gravito-electric propagation equation is 
given by 



Eab = -0^ah + CUrl{Hab) - \ {p^ + p^) Gab 



(32) 



where Hab is the "magnetic" part of Weyl tensor. The 
gravito-magnetic propagation equation is 

Hab = -OHab - curl{Eab) + ]^curl{Iiab) ■ (33) 

The system of equations is closed with the vorticity prop- 
agation equation which governs rotational perturbations: 



2 1 

Wa = -7:'diOa - -CUrl{Aa) 



(34) 



In what follows we will use the following linearized iden- 
tities for closed Einstein static background, for any scalar 
g, vector V"" and tensor Aab'- 

{y-g)-^V-{g), (35) 

(V"y^)- = V''(T>''), (36) 

(V^5)-=V2(5), (37) 



Va(V23) = v2(Va5)-^(Vag) 



(38) 



and 



,2.3, 3^ 



CUrl{cUrlAab) = -VMafc + — A„t + -V<aVM6>c . (39) 

a5 2 

Finally, to convert the system of partial differential equa- 
tions to ordinary differential equations, we use the stan- 
dard procedure of harmonic decomposition, employing 
the trace- free symmetric tensor eigenfunctions of the spa- 
tial Laplace-Beltrami operator defined by 



fc2 

v2q = — ^g, g = o. 



(40) 



We can then expand all first order quantities (both scalar 
and tensor) as 



x(i,x)-^x''^(t)g'=(x) 



(41) 



For the spatially closed models, the spectrum of eigen- 
values is discrete and given by [1^, 12^ 



fc2 ^n{n + 2) , 



(42) 



where the co- moving wavenumber n = 1,2,3,---. The 
mode n — Q corresponds to a change in the background. 
The first inhomogeneous mode n = 1 is a gauge mode 
which is ruled out by the constraint equations Goi — T^i , 
and reflects the freedom to change the 4-velocity of the 
fundamental observers. It follows that physical modes 
are constrained by fc^ > 8 [l^, |2^ . 

In the following sections we will extract the scalar, vec- 
tor and tensor components from the linearized pertur- 
bation equations using the standard procedure [28| and 
analyze the stability of the Einstein Static universe. 



A. Scalar Perturbations 

In order to derive the equations governing scalar per- 
turbations, we define variables describing the spatial 
Laplacians of density, expansion and spatial curvature 

A,„ = ^!ZV:, z = a^V^e, C - a^\I^R , (43) 

pm 

and the dimensionless Laplacians describing the inhomo- 
geneity in the Ricci scalar and its momentum 



n = a^V^R, ^^a^V^R 



(44) 



In what follows we obtain the evolution equations for 
the above scalar variables. Using the linearized identities 
together with (|28|) we get 



A. 



-{l + w)Z 



(45) 



The generalized Raychaudhuri equation pop and the 
trace equation ([9]) can be combined to eliminate the R 
term, giving 



(46) 



z= - 


R 

.3(l + w). 


Am 


4- 




— V^A^- 
iu -f 1 


'l + Sw 
2R 


v27^ 


Again from the linearized identities we get 


TZ^n, 


and using ^ we have 


2i?2 

"=9 


1 - 3w 

_{l + w){l+3w)_ 


A„ 


-1- 


R 
3 


3 
3 


w — 

w + 


1' 
1 


n 


+ \/^n 





(47) 



(48) 



Finally the constraint and the propagation equations for 
the curvature perturbation C are 



C_ 
7a 



\ ^^ 1 


A,„ + 


"l + Sw" 


.3(l + w) 


R 



v^7^ 



(49) 



and 



IV. VECTOR PERTURBATIONS 






'2R 
3 


Z + 


"l + Sw" 



V^3? 



(50) 



The system of equations p5]) - ((50)) are coupled partial 
differential equations. We use harmonic decomposition to 
convert them into a system of coupled first order ordinary 
differential equations. Following [321] we can then write 
them as a pair of coupled second order equations: 



2-e 



A 



(l+w)(6 + P+3ii;/c2) 



12 



TT 



(51) 



TV 



8(1 - 3u;) 



a4(l + 4w + 3u;2) 
6w - ^wk'^ - 2 - fc2 
a§(l + 3w) 



A* 



7^'= 



(52) 



where we have used R — 6/aQ. The above equations can 
be combined to a single fourth order equation for A^: 

^A^, - A'=(u;)^Af„ + B\w)^l, = , (53) 



where 



and 



A^(w) 



a2(l + 3w) 



(54) 



^(") = ^ 3^^I(rT3^io ^- ^''^ 

The general solution for equation ([55]) can be easily ob- 
tained as 

A^(t) = Cie-^-* + C2e+'^-* + C3e-^+* + C4e+^+* , (56) 

where 



/3± = W2yl'=(u;) ± 2./^'=(w7)2 - A,B^{w) 



(57) 



It is easy to see that the solution ((56)) is non-growing if 
and only if all the conditions 

1. A^iw) < 

2. B'^iw) >0 

3. A'^(w)2 -4B'=(u;) > 

are satisfied, and it follows that all three conditions are 
satisfied for all fc € [8, oo) if 



w > 



A + VE 
6 



(58) 



In other words, the Einstein Static universe is stable with 
respect to all scalar perturbations for Wc < w < 1. This 
range includes radiation and stiff matter, but not pres- 
sure free matter. It is also interesting that the ratio of 
Wrad and Wc is the golden ratio. 



The vector perturbations are sourced by the vorticity 
tensor tOab- At linear order, the vorticity propagation 
equation for the Einstein Static background is simply 



^ab 







(59) 



It follows that the Einstein Static model is neutrally sta- 
ble against vector perturbations for all equations of state, 
on all scales. 



V. TENSOR PERTURBATIONS 

Tensor perturbations [33] are characterized by the evo- 
lution of the transverse, trace-free part of the shear ten- 
sor ffab- At linear order the shear propagation equation 
decouples: 



o-afc = V^o-Qb + 7;ip'^ 



P 



R)o-ab 



(60) 



where we have only kept the tensor contributions. Sub- 
stituting the exact solution for the Einstein Static back- 
ground and decomposing into Fourier modes with co- 
moving index fc, we obtain 



C^ab 



P + 2 



CTafc = . 



(61) 



Solving this equation, it can easily be seen that for all 
equations of state, the Einstein Static model is neutrally 
stable against tensor perturbations on all scales. 



VI. CONCLUSIONS 

In this work, we have given a careful analysis of the 
stability of Einstein Static universe in fourth order the- 
ories of gravity with respect to generic inhomogeneous 
and anisotropic perturbations. 

First of all, we showed that the modified field equations 
only admit an Einstein Static solution for the very special 
form of gravitational Lagrangian: f{R) = a + bR'^, where 
the constants a, 6, c are fixed functions of equations of 
state parameter w and cosmological constant A if A ^ 0. 
If A = 0, a = and b is arbitrary. This is a surprising 
result, since one could naively expect that the additional 
degrees of freedom arising from a modified f{R) would 
make it easier to find solutions of the Einstein Static type. 
This result differs from some recent work 3q , where the 
stability of the Einstein Static solution was investigated 
for a number of types of f{R) theories which appear not 
to fall into the class discovered here. 

Our stability analysis shows that as in GR [1^, the 
Einstein static model is neutrally stable against vec- 
tor and tensor perturbations for all equations of state 
at all scales. Scalar perturbations are only stable on 



all scales if the matter fluid equation of state satisfies 
Cs > 'e~^ ~ 0.21. This result is remarkably similar to 
the GR case, where it was found that the Einstein Static 
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